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Abstract We searched for long period variation in V^-band, /^-band and RXTE X-ray 
light curves of the High Mass X-ray Binaries (HMXBs) LS 1698 / RX J1037.5-5647, HD 
110432 / IH 1249-637 and HD 161103 / RX J1744.7-2713 in an attempt to discover 
orbitally induced variation. Data were obtained primarily from the ASAS database and 
were supplemented by shorter term observations made with the 24- and 40-inch ANU 
telescopes and one of the robotic PROMPT telescopes. Fourier periodograms suggested 
the existence of long period variation in the V"-band light curves of all three HMXBs, 
however folding the data at those periods did not reveal convincing periodic variation. At 
this point we cannot rule out the existence of long term F-band variation for these three 
sources and hints of longer term variation may be seen in the higher precision PROMPT 
data. Long term V"-band observations, on the order of several years, taken at a frequency 
of at least once per week and with a precision of 0.01 mag, therefore still have a chance 
of revealing long term variation in these three HMXBs. 

Key words: accretion - stars: Be stars - stars: neutron stars - stars: individual: LS 1698, 
HD 1 10432, HD 161 103 - Xrays: binaries 



1 INTRODUCTION 

High Mass X-ray Binaries (HMXBs) consist of a massive star orbited by a compact object like a neutron 
star, a black hole or, as is proposed for 7 Cas analogues, a white dwarf. A major class of HMXBs, the 
BeXs, have a Be star as the primary massive star The orbital periods of many BeXs are known from 
their periodic Type I X-ray outbursts. Such outbursts are believed to be caused by the passage of a neu- 
tron sta r through the Be star's equato rial decretion disk when increased accretion onto the neutron star 
occurs (lOkazaki & Negueruela[|200l1) . These Type I X-ray outbursts can also be accompanied by optical 
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Table 1 Observation Log. 



Source 


Dates ^ 


Filters 


Telescope 


LS 1698 


21/12/2000 to 3/12/2009 


VJc 


ASAS 




24/8/2005 to 5/9/2005 


B,V,R,Ic 


ANU 40-in 




24/3/2006 to 1/4/2006 


B,V,R,Ic 


ANU 24-in 




17/7/2007 to 28/7/2007 


B,V,RJc 


ANU 40-in 




8/12/2008 to 12/2/2009 


VJc 


PROMPT-5 


HD 110432 


29/1 1/2000 to 29/7/2009 


VJc 


ASAS 


HD 161103 


31/1/2001 to 2/7/2009 


VJc 


ASAS 




4/8/2002 to 14/10/2002 


B,VJc 


ANU 24-in 




29/8/2005 to 5/9/2005 


B,V,RJc 


ANU 40-in 




23/3/2006 to 1/4/2006 


B,V,RJc 


ANU 24-in 




19/7/2006 to 19/7/2006 


B,V,RJc 


ANU 40-in 



^ day/month/year 



outbursts which can confirm the orbital period (ICoe &Edge[l2004l) . However, many BeXs including the 
three of interest here, are persistent X-ray sources, indicating a more uniform accretion rate and circular 
orbits with longer periods. The lack of regular X-ray outbursts means that the orbital periods for many 
of the persistent BeXs are unknown. In spite of the ex pected lack of corresponding optical outbursts , 
optical searches for orbital periodicity can be successful jCoe et al.L 120051 : ISchmidtke & Cowievil2006l) . 

Here we report on the search for long periods in the optical photometry of three persistent BeXs: LS 
1698, the opticalcounterparttoRXJ1037.5-5647:HD 110432, the optical counterpart to IH 1249-637 
and: HD 161103, the optical counterpart to RX J1744.7-2713. The source LS 1698 has been argued 
to be an analogue of X Per (^eig & Roche, 1999) while HD 1 10432 and HD 161 103 have been argued 
to be analogues of 7 Cas. X Per is a BeX with a neutron star and a long orbital period of 250.3 d 
(iDelgado-Martf et"aUl200l| ). The star 7 Cas has be en proposed to contain a white dwarf as its compact 
object, although single star models for 7 Cas exist (iLopes de Oliveria et al. L 120061 l2007l) . 

Th e main component of our data came from the All Sky Automated Survey (ASAS, IPojmanskH 
(Il998b ) supplemented with confirming data obtained with the 24 and 40 inch telescopes of the Australian 
National University (ANU) and with one of the 16 i nch telescopes of the Panchromatic Robotic Optical 
Monitoring and Polarimetry Telescopes (PROMPT iReichart et al.l (|2005|) ) array. 

This paper is organized as follows. In Section|2]we present the photometric observations, including 
RXTE X-ray photometry for LS 1698. Section |3]gives the results of formal searches for periods in the 
data using Fourier techniques. In Section |4] we give a brief review of what is currently known about 
the three HMXBs to motivate further investigation of these sources. Our conclusions are summarized in 
Section |5] 

2 OBSERVATIONS 

A total of four small telescopes/telescope systems were used to make optical observations of LS 1698, 
HD 1 10432 and HD 161 103. The observation log is given in Table [T] 

The bulk of the data were obtained with the ASAS-3 telescope system. ASAS-3 consists of two 
wide field 87-mm diameter f2.8 Minolta telephoto lenses, one equipped with a Johnson V filter, the 
other with the Johnson-Cousins Ic filter. The individual telescopes use Apogee AP-10 2048x2048 
pixel CCD cameras as detectors. ASAS-3 is an automated instrument and produces reduced photometry 
automatically. As part of the reduction process, the software rates the quality of the data. Only observa- 
tions rated as 'grade A' were used. The ASAS instrument is located at the Las Campanas Observatory 
in Chile. 

Two telescopes owned by the ANU and located at the Siding Spring Observatory in Australia were 
used. These telescopes had apertures of 24 and 40 inches. 

For the first observing run with the ANU 24-in telescope, from 4/8/2002 to 14/10/2002, for HD 
161103 a pulse-counting GaAs photo-multiplier tube photometer was used with a filter wheel con- 
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Comparisoii/Check stars 



Stai- 


RA 


Dec 


B 


V 


R 


Ic 


a 


10 37 33.4 


-56 47 19.4 


12.82 


12.45 


12.12 


11.71 


c 


10 37 47.17 


-56 51 29.2 


11.62 


11.17 


10.81 


10.41 


d 


10 37 23.9 


-56 51 11.4 


12.52 


12.05 


11.67 


11.28 


e 


10 36 52.05 


-56 48 20.9 


11.54 


11.18 


10.83 


10.41 



Fig. 1 Comparison/check stars in the field of LS 1698 as determined from data obtained 
with the ANU 24- and 40-in telescopes at Siding Spring Observatory. Errors are on the order 
of 0.01 to 0.05 mag. The target star is marked with two lines in the centre. Positions and 
magnitudes are listed for those stars whose magnitudes and colours make them most useful 
for differential photometry with the target. North is up, east to the left. 



trolled by a computer that also recorded photon counts in the B, V and Ic bands (iBesselL Il990l) . 
Transformations to the standard Cape U B VRI system were ma de on one night from 15 stars in the 
Regions E5, E7 and E9, using the values of iMenzies et al.l (Il989l) . On each night, the V , {B — V) and 
{V — I) magnitudes of LS 1698 were refered to the two standards E6/04 and E6/66, which are approxi- 
mately 18° from the source. The internal photometric errors of the observations were of the order 0.002 
to 0.004 mag. Referral to the fairly distant E6 Region stars adds an uncertainty of similar amount on 
average, depending on the air-mass and the quality of the night, due to the uncertainty in the differen- 
tial extinction corrections. The final differential photometric errors were therefore in the range 0.005 to 
0.007 mag. 

Subsequent observing runs with the ANU 24-in telescope, for LS 1698 and HD 161103, employed 
an imaging camera with a SITE 2048 x 2048 pixel LN2 cooled CCD and a manuall y turned filter whee l 
with B, V, R and Ic filters. The camera wa s controlled wit h the CICADA software ( Young et al.lll999l) . 
Standard stars in the Landolt SA104 field (iLandolll Il992h were observed for the purpose of defining 
comparison stars in the field of LS 1698. From the Landolt SA104 stars, extinction and filter band trans- 
formation coefficients were computed and combined with similar data from the ANU 40-in telescope 
(see below) to define the magnitudes of standard stars in the LS 1698 field (see Fig.[T]i- Standard star 
magnitudes in the field of HD 161103 were determined from data obtained at the Sonoita Research 



4 



G.E. Sarty et al. 





■ Nt 






n. . •■ ' 


• f 






°* • 






f 


ft- 








. E • • : 








" ■ ■ ■■■'■< 
■ .... ■*•■«•: 

h*- • • ■ 

















Comparisoii/Check stars 



Stai- 


RA 


Dec 


B 


V 


R 


Ic 


a 


17 44 41.9 


-27 14 31.5 


12.881 


12.082 


11.587 


11.053 


b 


17 44 44.0 


-27 14 22.5 


13.056 


12.590 


12.311 


12.002 


d 


17 44 52.1 


-27 16 10.5 


11.298 


10.731 


10.427 


10.081 


f 


17 44 12.9 


-27 11 11.5 


12.525 


11.936 


11.602 


11.228 


j 


17 44 28.4 


-27 18 57.4 


12.849 


12.222 


11.833 


11.392 


k 


17 44 40.2 


-27 17 33.7 


13.281 


12.648 


12.216 


11.817 


1 


17 44 47.3 


-27 14 05.5 


13.265 


12.716 


12.352 


11.983 



Fig. 2 Comparison/check stars in the field of HD 1 6 1 1 03 as determined from data obtained at 
the Sonoita Research Observatory (SRO) by Arne Henden of the AAVSO. Errors are on the 
order of 0.001 to 0.005 mag. The target star is marked with two lines in the centre. Positions 
and magnitudes are listed for those stars whose magnitudes and colours make them most 
useful for differential photometry with the target. North is up, east to the left. 



Observatory (SRO) and were provided by Arne Henden of the American Association of Variable Star 
Observers (AAVSO) (see Fig.|2]i. The errors on the magnitudes of the standard stars in the LS 1698 field 
were in the range 0.01 to 0.05 mag. The errors on the target star magnitudes, as determined by the inho- 
mogeneous differential photometry technique (iHonevcuttI 1 1 992h . were generally 0.005 mag (see Figs.|3] 
to |5]). Inhomogeneous differential photometry determines the differential magnitudes of all analysed 
stars at once, using a least squares approach, and references the differences relative to one selected star. 
For LS 1698, the reference star was star a as labelled in Fig. [l] For HD 161103, the reference star was 
star c as labelled in Fig.|2] The magnitudes of the reference stars were therefore added to the differential 
magnitudes as determined by the inhomogeneous photometry technique. 

The 40-in telescope was equipped with the Wide Field Imager (WFI, iBaade et al.Nl998h l and a 
filter wheel with B,V,R and Ic filters. The LN2 cooled camera and filter wheel were controlled with 
the CICADA software. Similar to the observing procedure described for the 24-in ANU telescope, stan- 
dard stars in the Landolt SA104 and SAl 13 fields were observed for the purpose of defining comparison 
stars in the source fields. As with the description given for the 24-in telescope above, extinction and filter 
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Fig. 3 (a) LS 1698 V band data, (b) LS 1698 Ic band data. A few oudier AS AS data points 
below V = 11.8 and Ic ~ 11-0 have been cut-off in the plots. 



band transformation coefficients were computed, comparison star magnitudes determined, and source 
magnitudes determined using the inhomogeneous photometry approach with similar errors. The differ- 
entially determined magnitudes of the target stars, for both the 24-in and 40-in telescope data, were 
not transformed to a standard filter band system since the precision of the transformation coefficients 
was deemed adequate only for determining comparison star magnitudes and not for use in differential 
photometry. As a check on the comparison star magnitude determinations, a comparison with magni- 
tudes published in the Simbad database was made where possible, with excellent agreement between 
the determined and published magnitudes. 

The fourth telescope used was the PROMPT-5 telescope located at the Las Campanas Observatory 
in Chile. The PROMPT telescope system is a robotic system consisting of five 16-in Richey-Chretien 
telescopes built by RC Optical Systems. Each telescope is equipped with a fast readout (< 1 s) Apogee 
Alta U47+ 1024x1024 CCD camera. The PROMPT system is under the control of Sky net, a prior- 
itized queue scheduling and control program running on a computer at UNC-Chapel Hill's Morehead 
Observatory. Skynet interacts with MySQL databases and commands dumb-by-design 'Terminator' pro- 
grams at each telescope. Images are automatically transferred back to a 1.1 terrabyte RAID 5 array at 
Morehead. 

Calibration of the CCD data for the ANU telescope and PROMPT data, involving bias subtraction, 
flat fielding from sky flats taken in morning and evening twilight and, for the PROMPT data, dark frame 
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Fig. 4 HD 110432 V and Ic band data. A few T^-band outliers below 7th magnitude have 
been cut-off in the plot. The data before HJD 2452300 were exposed with a different protocol 
than later data and non-linear CCD response for this bright source may be responsible for the 
apparently discrepant data points at the beginning of the light curve. 



subtraction, were performed with the IRApQ software and the ccdproc task. Instrumental magnitudes 
for the ANU t elescope and PROMPT data were determined with the point-spread-function approach 
(IStetsonL[l987l) using the psf and allsky tasks of the IRAF software. The instrumental magnitudes were 
then subjected to the inhomogeneous differential photometry approach as previously described. 

Data taken for the three sources in the V and Ic bands are shown in Figs. [3]to|5] Data from the 
B and R bands were not of sufficient temporal density to be useful for period determination; they are 
summarized in Table|2]for completeness sake. 

X-ray data forRX J 1037.5 —5647, for which LS 1698 is the optical counterpart, were available from 
the Rossi X-Ray Timing Explorer (RXTE) satellite. Definitive data (as opposed to quick-look data) from 
the All Sky Monitor (ASM) instrument were downloaded from the RXTE Guest Observer Facility for the 
purposes of finding correlations with the optical wavelength variations. The RXTE ASM is composed of 
three Scanning Shadow Cameras (SSCs) that p erform sets of 90 s p ointed observations (dwells) covering 
about 80 per cent of the sky every ~90 min dLevine et al l fl996l) . The data used were from the dwell 
by dwell compilation. Each raw data point represents the fitted source flux from one 90 s dwell. Data 
from all three SSCs were used and represent nominal 2-10 keV rates in ASM counts s~^. Nominally, 
the Crab nebula flux is about 75 ASM counts s^^ (when the source is at the centre of an SSC field of 
view and all eight anodes are operational). We retained only the '3cr' detections based on the reported 
variance for each data point. The ASM data for LS 1698 / RX J1037.5-5647 are shown in Fig.|6l 



' IRAF is the Image Reduction and Analysis Facility, a general puipose software system for the reduction and analysis 
of astronomical data. IRAF is written and supported by the IRAF programming group at the National Optical Astronomy 
Observatories (NOAO) in Tucson, Arizona. NOAO is operated by the Association of Universities for Research in Astronomy 
(AURA), Inc. under cooperative agreement with the National Science Foundation 
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Fig. 5 (a) HD 161 103 F band data, (b) HD 161 103 /c band data. A few AS AS outliers below 
V = 8.8 have been cut-off in the plot. 



Table 2 Summai-y of B and R band data from LS 1698 and HD 161 103. 



Source 


band 


Time span^ 


Mean 
mag 


Std. 
Dev. 


No. 
Obs. 


LS 1698 


B 


24/8/2005 to 5/9/2005 


12.004 


0.015 


8 






24/3/2006 to 1/4/2006 


11.947 


0.048 


6 






17/7/2007 to 28/7/2007 


12.043 


0.012 


26 




R 


24/8/2005 to 5/9/2005 


10.983 


0.010 


8 






24/3/2006 to 1/4/2006 


10.976 


0.053 


10 






17/7/2007 to 28/7/2007 


11.037 


0.009 


29 


HD 161103 


B 


4/8/2002 to 14/10/2002 


8.891 


0.017 


29 






29/8/2005 to 5/9/2005 


8.898 


0.135 


10 






23/3/2006 to 1/4/2006 


9.223 


0.027 


3 






19/7/2006 to 19/7/2006 


8.904 


0.006 


11 




R 


29/8/2005 to 5/9/2005 


8.041 


0.017 


11 






23/3/2006 to 1/4/2006 


8.081 


0.037 


9 






19/7/2006 to 19/7/2006 


8.058 


0.010 


11 
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LS 1 698 - RXTE data 




2000 3000 4000 

HJD - 2450000 



Fig.6 2-10 keV RXTE light curve for LS 1698 / RX J1037.5-5647. The red Une indicates 
the time-span of of the optical observations shown in Fig.|3] 



3 RESULTS 

Fouri er analysis of the photometric data was performed using the software PERIOD04 jLenz & Bregeri 
12005^. Analyses were performed with either the data weighted by the inverse of the errors or with 
equally weighted data. Unless otherwise reported (e.g. for the RXTE data) analyses with data weighted 
by the inverse of the errors did not produce any sensible results. The largest peak was deemed significant 
if the signal to noise ratio (SNR), as computed by PERIOD04, was greater than 4. When such significant 
peaks were found the data were folded at the indicated period and two Kolmogorov-Smirnov (K-S) tests 
were performed. The first K-S test tested the null hypothesis that the light curve had a constant value. 
The second K-S test tested the null hypothesis that the light curve data had a Gaussian distribution with 
the mean and standard deviation of the data. For each K-S test, we report the probability, p, that the 
data represent the null hypothesis and the distance between cumulative distributions, d, upon which 
the probabi l ity is based. All K-S statistics were computed using the Numerical Recipes algorithms 
jPress et al.Lll986l) . 

For the source LS 1698, period searches were done for the entire V- and /c-band datasets. It was 
not necessary to trim the obvious outliers (see Fig. |3} from the data prior to analysis since analysis with 
arbitrarily defined outliers removed did not change the results. The search results are given in Table |3] 
and the resulting periodograms are shown in Fig. |7] Note the discrepancy between the periods found 
in the V- and /c-bands. A visual inspection of Fig. |7] (look at the window function shape) and Fig. [3] 
suggests that the higher frequency found for the /p-band data is likely an artifact of the sampling gaps 
and that the variation of the /p-band data roughly correlates with the variation of the l^-band data. For 
example, a minimum for both light curves occurs around HJD 2453800. We therefore regarded the 645 
d period as a viable period for LS 1698. Folding the data at the 645 d period did not reveal a convincing 
periodic variation, see Fig.Qc). The first K-S test did not show any evidence of variation, d — 0.002, 
p = 0.99. The second K-S test indicated that the variation of the data was significantly different from 
Gaussian noise, d = 0.170, p = 1 x 10~^^. The more densely sampled PROMPT-5 data, which show a 
hint of periodic variation upon visual inspection, were subject to a separate search for higher frequencies 
but no significant periods were found. A period search of the unweighted RTXE data for periods yielded 
a period of 184 d at a m arginal SNR of 4.2. The period 184 d is almost exactly a half year and is likely 
an instrumental artifact (IWen et al.l 120061; ISartv et all l2009h . A search with weighted data showed no 
significant periods. 
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Fig. 7 (a) Periodogram for the LS 1698 V band data. The peak frequency corresponds to a 
frequency of 0.00155 cycles d^^, which is a period of 645 d. (b) Periodogram for the LS 1698 
Ic band data. The peak frequency corresponds to a frequency of 0.001 18 cycles d^^, which 
is a period of 847 d. (c) V band data folded at the 645 d period. No consistent variation at 
the 645 d period is apparent, however note that the higher precision PROMPT data hint at a 
possible shorter period variation in the light curve. 
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Fig. 8 Periodogram for the HD 110432 V band data. The peak frequency corresponds to a 
frequency of 0.00054 cycles d~^, which is a period of 1852 d. 



Table 3 Period search results. 



Source 


Band 


Period 
(d) 


Amplitude 
(mmag) 


SNR 


LS 1698 


V 


645 


35 


7 




Ic 


847 


155 


18 


HD 110432 


V 


1852 


134 


7 


HD 161103 


V 


1075 


37 


7 




Ic 


532 


52 


7 



For the source HD 1 10432 a formal period search was done for the T^-band dataset. Outliers did not 
need to be removed from the data prior to analysis. The search results are given in Tableland the result- 
ing periodograms are shown in Fig. |8] The most significant period of 1852 d essentially corresponds to 
the length of the data set, so continued observations for at least another 5 years are required to confirm 
a period of that length. 

For the source HD 161 103, period searches were done on the entire V- and /p-band datasets, again 
without trimming the small number of apparent outliers from the datasets. The search results are given 
in Table|3]and the resulting periodograms are shown in Fig. |9] Folding the data at the 1075 d period did 
not reveal a convincing periodic variation. The first K-S test did not show any evidence of variation, d = 
0.002, p — 0.99. The second K-S test indicated that the variation of the data was significantly different 
from Gaussian noise, d = 0.130, p = 1 X 10"". Once again, there was a discrepancy between the 
periods found for the V- and /p-band data. The higher frequency period found in the /c-band data may 
again be an artifact of the gaps in the sampling. A visual inspection of Fig. |9]reveals a possible negative 
correlation between the V- and the /p-band light curves. We will comment on this anti-correlation in 
Section m 





Fig. 9 (a) Periodogram for the HD 161 103 V band data. The peak frequency corresponds to 
a frequency of 0.00093 cycles d~^, which is a period of 1075 d. (b) Periodogram for the HD 
161 103 Ic band data. The peak frequency corresponds to a frequency of 0.00188 cycles d~^, 
which is a period of 532 d. 
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4 DISCUSSION 

As a part of motivating further long-term optical observations of LS 1698, HD 1 10432 and HD 161 103, 
we give a brief background of what is known about these three sources. 

4.1 Background for the three systems 

4.1.1 LS 1698/RXJ1037.5-5647 

The source 4U 10 36-56 (believed to be the s ame as LS 1698 / RX J1037.5 -5647 according to 
iMotch et al] (1 19971) ) was first observed by Uhuru dForman et al.Lll978l) and OSO 7 jMarkert et al.l 1 19791) 
with a mean Uhuru flux of 1.0 x 10^^" erg cm^^ s~^ (2-10 keV). A flare appeared in November 1974 , 
with a flux of 2.4 x 10^^° erg cm^^ s^^ (2-10 keV) in observations made bv Ane/ V (IWarwicklll98ll) . 
/?(95Ar detected a flux about 1 times less than the Uhuru flux between 1970 and 1976 with later pointed 
RXTE observations by lMotch e t al. ( 1997) being about 20 times fainter than the ROSAT 1970 - 1976 val- 
ues. The optical counterpart, LS 1698, of RX J1037.5-5647 was identified bv lMotch etal] d 19971) on 
the basis that i t was within the 95% confidence region of the ROSAT error circle from the first ROSAT 
afl-sky survey dVogeslI 19921: Motch et al.lll99ll) and that Lx/Lboi > 10^^ which was consistent with 



being an X-ray binary. IMotch et al.l (119971) further obtained blue and red spectra with the EFOSC2 in- 



strument on the ESO-MPI 2.2-m telescope that indicate a BO V-IIIe classification with a reddening of 
E{B — V) = 0.75±0.25, meaning a probable distance of ~5 kpc as opposed to an earli er estimate of 
18 kpc. An X-ray pulse period of 860±2 s has been measured bv iReig & Rochd d 19991) for LS 1698 
/ RX J1037.5-5647 in RXTE Proportional Counter Array (PCA) observations with X-ray flux levels 
consistent w ith the earlier Uhuru a nd Ane/ V levels. Assuming a distance of 5 kpc, the X-ray flux levels 
measured bv lReig & Rochd d 19991) translate to an X-ray luminosity of ^--^4.5 x 10'^^ erg s~^ (3-30 keV). 
The existence of X-ray pulses and a 10^^ e r g s~^ flux level indicates the presence of a spinning, accret- 
ing magnetic neutron staj^ and the ICorbetlJl986l) spin-orbit period relation for BeXs suggests that the 



orbital period for LS 1698/RX J1037. 5-5647 is on the order of 200 d. Our data indicate that the orbital 
period is possibly 645 d. The X-ray variability appears to be around a factor of 10 with the observed 
flares possibly being caused by sporadic ejections of Be star disk material; with such unsteady accretion 
we might expect that the neutron sta r spin has not reached an equilibrium and therefore the system will 
not fit the Corbet spin-orbit relation (lReig&Roche[|l999b . The faster spin may be due to an equilibrium 
with the accretion of the denser clumpy material instead of with the background density of material at 
the the neutron star's orbit. 



4.1.2 HD 110432/ IH 1249-637 

HD 110432 (= BZ Cru = HR 4830), a membe r of th e 60 Myr old open cluster NGC 4609 be- 
hind the southern Coalsack dFeinstein & Marracot fl979l) . was id entified as the optic al counterpart to 



IH 1249-637 bv iTuohv et al.l dl988l) . Optical observations bv iDachs etal] d 19891) with the 61-cm 
University of Bochum Cassegrain telescope at La Silla, Chile supported an MK classification of B 1 
Ille. Model fits to their optical spectra indicate Toff — 22500 K, log g ~ 3.5, i?* = 8.6 R., a projected 
rotation velocity Vyot sin i ~ 300 km s~^, an interstellar extinction of E^^{B — V)~ 0.34 and a circum- 
stellar extinction of E'^^ {B — V) =0.16. UV observations with the lUE space telescope bv lCodina et al.l 
(fr984i) indicate a rotational Wjot siiii for the primary of 360 km s ^, T^s = 25000 K, \ogg = 3.5, a 
reddening of E{B — V)= 0.40 and spectral classific ation of B0.5 Ille with an absolute magn itude of 
Mv = —4.3 in general agreement with the results of iDachs et al.l (Il989l) . ICodina et al.l d 19841) quote a 
distance of 430 ± 60 pc which is in good agreement with HIPPARCOS parallax measurements that give 



^ A neutron star radius and mass, R = 10^ cm and M = 1 .4 Mq , ai'e implied by Lx = rjGihM / R for reasonable accretion 
rates, rh ~ 10~^^ Mq yr~^ and ~ 10'^^ erg. Here G is the universal gravitational constant and rj ~ 0.5 is the efficiency 
with which gravitational potential energy is converted to light. Pulsations would be caused by the magnetic funnelling of accreting 
material onto magnetic pole hot spots on the neutron star that rotate in and out of view. 
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the distance to HD 1 10432 as 300^40 pc (IChevalier & Ilovaisky[ll998h . The magni tude and pos ition of 
the Balmer discontinuit y implies T^s = 225 10 K, \ogg = 3.9 and M* = 9.6Mq (IZorec et all [2005jl). 
The UV observations of ICodina et al.l (1 19841) further show broad photospheric absorption lines, narrow 
emission lines of N V, C IV and Si IV with radial velocities of —1350 km s^^ and asymmetric profiles, 
for C IV and Si IV, that suggest an expanding envelope corresponding to a mass-loss of 3 x 10~^ Mq 
yr~^. The rotational velo city of v^ot sini = 300 km s^^, if matched with the 1.77 d photometric period 
found by |Barreral(ll99lb and the ra d ius of a BO. 5 III star, 1 1.55 Rq, indicates a primary spin axis in- 



clination of i = 58° ( lBalonaill995b . [Smith & Balonal (l2006l) did not obser ve the 1.77 d period in their 



optical photometry but report a possible 130 d period. ICodina etalJ(ll984 also note a strong polariza- 
tion of the light from HD 1 10432 in observations made with a photopolarimeter attached to the 1.6-m 
telescope of the Brazilian Astrophysical Obsei-vatory (OAB) in Brasopolis, however they attribute most 
of the polarizatio n to the star's pro ximity to the southern Coalsack dark nebula which is between the 
star and the Earth jCrawford[[l99TI) . Radial velocities were obs erved to vary between 2 7 and 70 km s 



by iBuscombd ( Il962 ). li kely based on moving e mission lines: ISmith & Balonal ( l2006h adopt the latest 
value of 6 km s~^ from iThackerav et al.l (Il973l) . which is confirmed by their own measurements, for 
their work. 

Variable sub-features, wind absor ptions over emission lines, have been obse rved the UV (C IV 
A1548, Si IV A1403 and N V A1238) dCodina et aUll984l:ISmith & Balonall2006l) . More dramatic mi- 
grating subfeatures (msf's) have been observed in the optical Hel A5876 and A6678 emission lines 
with accelerations on the order of 100 km s~^ hr~ -1 dSmith & Balonal |2006|) . The msf's have been 
hypothesized to be due to strong magnetic connections between the Be star and its decretion disk 



dStnith & Balona , 20061) and associated magnetic recon nection events may be the source of X-ray emis 



sion ( Smith & Robinsoiil 1 1999 1: ^Robin son et al.L |2002|) . The optical emission lines of Fe II and He I 



display double-lobed profiles typical of a disk, the Be star's decretion disk, seen nearly edge-on. Using 
spectral synthesis techniques. ISmith & Balonal (|2006|) find the temperature and density of the disk to be 
roughly 9800 K and 3 x 10^^ cm~^ respectively. They also find that the projected disk covers a large 100 
stellar areas out to a distance of 1 AU (~30 R*). The inferred disk volume and mass are approximately 
10^^ cm'^ and 10^^ Mp). The strongest abso rption wings in the optical and UV lines extend out to at 
leas t ±1000 km s~^ dSmith & Balonall2006l) which may be due to circumstellar material. 

iTorreion & OrrI (l200lb report the detection of a 14 ks X-ray pulsation period from BeppoSAX ob- 
servations that they attribute to the partial occultation of a hot spot on the compact object, likely a white 
dwarf, as it spins. Later observations with XMM-Newton do not show a coherent 14 ks period but time - 
series of the 2-12keV/0.6-2keV hardness ratio do show a 14 ks period (iLopes de Oliveria et al I I2OO7I) . 
An equilibium between accretion angular momentum and spin for a white dwarf with a magnetic field 
of 10^ G and a spin period of 14 ks, implies an orbital period of several hundred days ([Apparao, 1994): 
our observed F-band period of 1852 d is longer by roughly a factor of two. However, based on the ob- 
servations described in the previous two paragraphs, HD 1 10432 appears to have a dense and complex 
circumstellar environment with a large dusty disk that interacts with the primary's magnetic fields. So 
the densi t y of a ccreting matter at the white dwarf orbit may be higher than that implicitly assumed by 



lApparaol ( 19941) and thus consistent with a 1852 d orbital period and an equilibrium spin period of 14 
ks. 

The flux observed by BeppoSAX was L2_iokcV = 3.4 x 10^^ erg s^^ assuming the HIPPARCOS 
derived distance of 300 pc. Past observations of IH 1249-637 with HEAO 1 and ROSAT (it is too faint 
for the RXTE ASM) indicate a persistent low luminosity source consistent with a wind-accreting white 
dwarl0 jTorreion & Onll200ltlmters[ll989l) . IH 1249-637 shows an emission feature at 6.8 keV with 
an equivalent width of 600 to 700 eV which is too large to b e associated with an ac creting neutron 
star but is similar to that found in some cataclysmic variables dTorreion & Orri 12001 1) , which contain 
white dwarfs. Observations with XMM-Newton show the 6.8 keV feature to be an Fe Ka complex of 
three lines (at 6.4, 6.7 and 6.97 keV) plus there is a suspected Fe XXVI Ly/3 emission at 8.2 keV 



' An accretion rate of m = 4 X 10 Mq yr ^ is implied by Lx = 3.4 X 10''^ erg s ^, rj = 0.5 and a white dwarf with 
10. 



10^ cm and M = 1.0 Mr. 
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dLopes de Oliveria et alll2007h . lTbn-ej6n & OrdfcOOll) note that E{B - V) ^ 0.40 implies an absorp- 
tion of Nh ~ 0.3 X 10^^ cm~^ but Nh ~ (1.1-1.4) x 10^^ cm~^ is deduced from the X-ray spectral 
model fits which imples a higher density of circumstellar material around the X-ray source than around 
the optical source. The existence of denser material around an accreting compact object may, again, be 
responsible for a relatively fa ster spin period of 14.5 ks at a relatively slower 1852 d orbital period. 
iLopes de Oliveria et al.l ( |2007|) find that a three thermal plasma model best fits their XMM-Newton X- 
ray spectra; the three plasmas have kT temperatures of [0.2-0.7 : 3-6 : 16-37] in a [few% : 15-25% 
: 70-80%1 distribution. The so ft component is consistent with a normal B star shocked wind emission 
(iLopes de Oliveria et al.l I2OO6I) . The hard component qualifies it, along with the low X-ray luminosity 
and the Fe Ka complex, as a 7 Cas analogue. The three components of the X-ray spectra again attest to 
the complex circumstellar environment of HD 1 10432. No coh erent 14.5 ks variations were fo und, but 
the XMM-Newton 0.6—2 and 2-12 keV light curves examined bv lLopes de Oliveria et aL I (120071) showed 
'flaring' activity characterized by a '^/^"^^ power spectra. The brightness of the flares was not corre- 
lated with their X-ray hardness. 



4.1.3 HD 161103 /KXJ1744.7-2713 



HP 161103 was id entified as the optical counterpart of the BeX candidate RX J1744.7-2713 by 
iMotch et alJ (Il997h . The identification was made on the basis of optical spectroscopic observations 
with the ESO-MPI 2.2-m telescope, which showed Balmer emission lines, and pointed observations 
with ROSAT io follow up observations from the RO SAT Ga\&cl\c Plane Survey. HD 161103 is a B0.5 
V-IIIe star with a hard (7 Cas-li ke) X-ray emission with an unabsorbed luminosity of ^10^^ erg s~ ^ 
and ix/iboi of about 4 x IQ-^ jMotch et al.Lll997l:ISteele et all 1 1999t ILopes de Oliveria et al.l l2006l) 
Similar to HD 1 10432, its X-ray spectra as observed with XMM-Newton dLopes de Oliveria et al.ll2006l ) 
could be equally well-fit with a optically thin hot plasma, with T ^10^ K and solar abundances (XSPEC 
mekal model), or by a power law with photon index F ~ 1.5-1. 8. The X-ray spectrum also shows an Fe 
Ka complex with a neutral fluorescent line at 6.4 keV, an He-like Fe XXV emission at 6.7 keV and an 
H-like Fe XXVI emission at 6.97 keV. Such Fe Ka complexes are not seen in BeXs where there is direct 
accretion onto the surface of a neutron star and the 6.4 keV flouorescence indicates a reprocessing of ra- 
diation in cooler circ umstellar matter. The H- and He -like Fe lines indicate a hot plasma compatible with 
the overall mekal fit dLopes de OUveria et al.l I2OO6I) . A two component model with hot and cold com- 

S ionen ts also fits the spectrum well with the hot component being less absorbed. iLopes de 01iveria "et al.l 
20061) observed a coherent sinusoidal oscillation with a period of 3245±350 s over most of a 17.6 ks 
observation in the 0. 5-12 keV band t hat may be related to the spin of a compact object. If so then the 
spin-orbit relation of lApparaol (1 19941 ) for Be/white dwarf systems would imply an orbital period of a 
few hundred days and, again, we would explain the longer period of 1075 d implied by our data as 
being due to a denser circumstellar environment than implicitly assumed by Apparao. A white dwarf 
magnetic field less than the 10^ G assumed by Apparao would also lead to a longer orbital period for 
a given white dwarf spin rate. In assessing any proposed spin-orbit relations for both HD 110432 and 
HD 161 103 it should be remembered that the X-ray pulsations were observed only once for each object; 
both spin periods still need further observation to confirm their existence. 

The BO. 5 classification alo ng with a luminosity class of V or III i mply a distance to HD 161103 
of 1.1 or 2.0 kpc respectively dLopes de Oliveria et al.l l2006t IWegneii [1994; Humphrevs & McElrovl 
Il984l) . The optical spectra of HD 161103 show a large equivalent width Ha line, -31 to -34 A, with 
weak V/R asymmetry along with the Pasc hen fine series in emission in the IR spectra that indicates 
the presence of a large stable decretion disk dLopes de Oliveria et al.ll2006l) . The l arge disk again p oints 
towards larger possible accretion onto a white dwarf than impUe d by the relation of lApparaol d 1 9941) . The 
O IA8446 line, excited by Ly/3 photons ( Andrillat et al.l Il988l) . is also seen in emission. An EW(Ha)- 
Porb relation discovered by Reig et all dl997 ) implies Porb > 100 d. A number of Fe emission lines 
from the circum stellar environment are visible while the He I A43 88 photospheric line varies by less 
than 13 km s"^ dLopes de Oliveria et al.l l2006l). ISteele et al.1 dl999b find v^ot sini = 244 ± 33 km s^^ 
from the He I lines. iLopes de OUveria et al.1 d2006b note that the strength of the N II A3995, 4044 lines 
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indicate that HD 161103 is moderately Nitrogen-rich, a quality seen in other HMXB primaries that 
could be the result of mass transfer in binary evolution - a fact that argues for the binary nature of this 
7 Cas analogue. 

4.2 Further searches and possible causes for orbitally induced optical variations 

Based on previous observations of the these, and other HMXBs, all three of the HMXBs we observed 
are expected to have long orbital periods, if they are, in fact, binary systems. Although we are unable 
to confirm any periods in the optical data we have so far, further long term observations, at a precision 
of 0.01 mag or better (e.g. at the the precision of our PROMPT data), may yet reveal periodic variation. 
If such variation is driven by orbital motion, it would raise the question of what physical mechanism 
is behind the the light curve variation. Ordinarily, two mechanisms are responsible for orbital period 
light curve variation: eclipses and ellipsoidal variation. Variation due to eclipses would be one or two 
localized dips in the brightness of the star, due to the primary and secondary eclipses, and not the smooth 
variation seen in our data. Ellipsoidal variation, caused by the primary being distorted into an ellipsoid- 
like shape by gravitational tidal forces, generally requires a short orbital period or a highly eccentric one 
with a small periastron distance. The orbital periods for the three sources studied here here are likely 
not short and a close periastron would result in regular X-ray outbursts which are not seen. Optical 
modulation can occur in synchrony with X-ray modulation with both bands brightening, possibly with 
the optical brightening occurring first, in response to i ncreased mas s flow through an accretion disk as 
the compact object passes perihelion dSchmidtke & Cowley. 2006; Kriss etaU[l983l) . Where we have 
X-ray data, for LS 1698, there is no correlation of the X-ray output with the optical light curve, ruling out 
simple accretion disk models in that case. The X-ray luminosities of the other two sources HD 1 10432 
and HD 161 103 are too low to produce a long-term record in RXTE data. So it may yet be true that the 
low luminosity X-ray light curves of HD 1 10432 and HD 161 103 coiTelate with their optical light curve. 

In the absence of the modulation of the optical light curve by direct accretion onto the compact 
object, a mechanism other than eclipses and ellipsoidal variation would have to be responsible for the 
variation seen in our data for LS 1698, HD 110432 and HD 161103. For example, decretion disk ac- 
tivity or the interactio n of the decretion disk w ith a white dwarf as causes for "type 1" Be star activity 
could be responsible dMennickent et al.l 1200 ih . However such type 1 activity generally takes the form 
of an optical outburst in contrast with the slower variation seen in our data. We therefore investigated 
two, speculative, geometric models for orbitally induced optical variation here without delving into the 
detailed physical mechanisms behind the models. With the geometry narrowed down by at fit of these 
models to a future data set we may then speculate about the required physical mechanisms. The two 
models we investigated are the 'emitting shock model' and the 'hot surface model'. 

In the emitting shock model it was assumed that a shock preceding the neutron star (or a possible 
white dwarf for the cases of HD 110432 and HD 161103) emits light by some unspecified physical 
process. Colliding winds from the compact star and the primary star may give rise to such a shock or the 
shock might be at the boundary of the compact object's magnetosphere. Generally such shocks would 
be expect ed to emit hig h energy 7 or X radiation (see e.g. models for the HMXBs LS 5039 and LSI 
H-61°303 (ICeruttil[2008h ) but here we allow for the possibility of optical band emission. For simplicity, 
the emitting shock was modelled as a half-lit surface of a sphere aligned with the orbital velocity vector 
of the neutron star. The half-lit surface would thus present a well-defined phase to the observer with the 
maximum of the light curve corresponding to the full phase, or at least maximum visible phase, that 
occurs when the neutron star is approaching the Earth. 

In the hot surface model it was assumed that the half of the primary star facing the neutron star 
is brighter than the other half. In this model the light curve maximum occurs when the neutron star 
is between the primary star and the Earth when the hot surface presents a full phase, or a maximum 
visible phase, to the Earth. To re-emphasize, the model is purely geometric; there is no modelling of 
brightness as a function of the inverse square distance between the component stars for example. A 
simil ar hot surface model, with an X-ray shadowing disk , has been proposed for the X-ray binary Her 
X-I (iGerend & BovntonLll976trHowarth & WilsonLll983h . 
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Both models produced light curves that visually match the variation seen in the unfolded light curves 
of Figs. [3]to|5] However, until convincing values for the orbital periods are found, it makes little sense 
to formally fit the two geometric models to the data. 

The variation in the V^-band magnitude of the lightcurves implies a variation in the energy outputs 
of the systems that we may estimate given the distances to the sources. Specifically we may estimate the 
luminosity of the system i* at maximum and minimum brightness via 



and take the difference to estimate the luminosity change. Here Lq, cIq and V© are the solar luminosity, 
distance and apparent l^-band magnitudes respectively; L^, and K are the same parameters for the 
system and = 3.2E{B — V) is the interstellar absorption. Using quantities found by others, as 
reported in the background sections above, we find a change in energy output of 1.3 x 10"^^ erg s^^ for 
LS 1698, 9.6 X 10^^ erg s"i forHD 110432 and 1.1 x 10^^ erg s"! forHD 161103. These values very 
nearly match the observed X-ray luminosity for LS 1698 and are brighter than the X-ray luminosity 
of the other two systems. This suggests that the X-ray output of the system is reprocessed to optical 
wavelengths when the sources are optically bright, if an argument can be made that nearly all of the 
X-ray output is intercepted by the reprocessing medium. In the case of HD 110432 and HD 161103, 
additional, or alternate, sources of energy are also implied. This difference in energy source between 
LS 1698 and the pair HD 1 10432 and HD 161 103 further supports the notion that the character of the 
systems are different, recalling that LS 1698 likely contains a neutron star and the other two systems 
likely contain a white dwarfs. 

In order for the X-ray energy to be reprocessed to optical wavelengths, an optically thick repro- 
cessing media is necessary. An emitting shock is likely to be optically thin and, if that is the case, it 
would rule out such a shock as a reprocessing source of the systems' variable optical output. The hot 
surface, if on the primary, would have to be on the cooler surface of the star. Another possibility for a 
reprocessing site could be protrusions of a cold disk. It may be possible, for example, that the geometry 
of an illuminated cold disk edge is similar to that of our hot surface model. 

Simpler explanations for the observed optical variation could be related to the expected dusty cir- 
cumstellar environments of BeXs. The primary Be stars likely have equatorial dusty decretion disks, 
as indicated in their optical spectroscopy (see the backgrounds sections above). The interaction of the 
compact object with that dusty circumstellar environment may produce the observed optical variations 
as denser regions near the neutron star (or white dwarf) obscure the light from the primary. This obscur- 
ing dust model has some support from the V- and /p-band light curves of HD 161 103 (see Fig. [S]); as 
the y-band light curve gets brighter, the /^-band light curve appears to dim. An obvious explanation is 
that as the dust clears, the source gets brighter in the V-band and as the dust obscures the primary we 
see that dust as increased /c-band brightness. 

Clearly, more optical band data, of a higher precision than those of the ASAS data analysed here, 
are required before any of the models outlined here can be considered further 

5 CONCLUSION 

Long term variation is evident in V-band light curve data for three HMXBs. Using distances determined 
by others, the variation in optical output from the three HMXB systems implies a variation in luminosity 
of ~^10"^^ erg s^^ in all three cases. This variation in luminosity matches the X-ray luminosity of LS 
1698 which is suspected of containing a neutron star The X-ray luminosity reported by others for HD 
1 10432 and HD 161 103, which are suspected of containing accreting white dwarfs, is smaller at ^10'^^ 
erg s~^. The light curve variations are roughly consistent with a hot surface source on the hemisphere of 
the primary facing a compact accreting object or with a hot shock-like structure preceding the compact 
object in its orbit. The variation may also, or instead, be due to obscuration by circumstellar dust and 
gas as denser regions follow the orbit of the compact object. The hypothesis of variation due to dust 




(1) 
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is supported by the observation of apparently anti-correlated V- and /c-band light curves in the case 
of HD 161103. A V- Jc-band anti-correlation is not apparent for the other two sources. The lack of 
periodicity in the 2-10 keV X-ray light curve for LS 1698 rules out a simple accretion disk brightening 
model for that source. 

ACKNOWLEDGMENTS 

The Digitized Sky Surveys used to produce Figs.[T]and|2]were produced at the Space Telescope Science 
Institute under U.S. Government grant NAG W-2166. The images of these surveys are based on photo- 
graphic data obtained using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt 
Telescope. The plates were processed into the present compressed digital form with the permission 
of these institutions. The Beowulf cluster used to investigate the geometrical models was constructed 
from recycled computers donated by the Arts and Science Computer Laboratory of the University of 
Saskatchewan. Thanks to Ame Henden of the AAVSO for the SRO data used to define standard stars 
in the field of HD 161103. Thanks to the ANU technical support staff and the wonderful cooks at the 
ANU lodge for their support and hospitality during the observing runs with the 24- and 40-in telescopes 
at Siding Spring Observatory. 

References 

Andrillat Y., Jaschek M. & Jaschek C, 1988, AASS, 72, 129 
Apparao, K.M.V., 1994, A&A, 291, 775 

Baade, D., Meisenheimer, K., Iwert, O., et al., 1998, Messenger, 93, 13 
Balona, L.A., 1995, MNRAS, 277, 1547 

Barrera, L.H., Mennickent, R.E. & Vogt, N., 1991, AAS, 185, 79 

Bessel, M.S., 1990, PASP, 102, 1181 

Buscombe, W., 1962, MNRAS, 124, 189 

Cerutti, B., Dubus, G. & Henri, G., 2008, A&A, 488, 37 

Chevalier, C. & Ilovaisky, S.A., 1998, A&A, 330, 201 

Codina, S.J., de Freitas Pacheco, J.A., Lopes, D.F., & Gilra, D., 1984, A&AS, 57, 239 

Coe, M. & Edge, W., 2004, MNRAS, 350, 756 

Coe, M., Edge, W., Galache, J. & Mcbride, V., MNRAS, 356, 502 

Corbet, H.D.C., 1986, MNRAS, 220, 1047 

Crawford, I., 1991, A&A, 246, 210 

Dachs, J., Poetzel, R. & Kaiser, D., 1989, A&AS, 78, 487 

Delgado-Marti, H., Levine, A.M., Pfahl, E. & Rappaport, S.A., 2001, ApJ, 546, 455 
Feinstein, A. & Marraco, H.G., 1979, AJ, 84, 1713 

Forman, W., Jones, C, Cominsky, L., Julien, P., Mumy, S., Peters, G., Tananbaum, H. & Giacconi, R., 1978, ApJS, 
38, 357 

Gerend, D. & Boynton RE., 1976, ApJ, 209, 562 
Honeycutt, R.K., 1992, PASR 104, 435 
Howarth, I.D. & Wilson, B., 1983, MNRAS, 202, 347 
Humphreys, R.M. & McElroy, D.B., 1984, ApJ, 284, 565 

Kriss, G.A., Cominsky, L., Remillard, R.A., Williams, G. & Thorstensen, J.R., 1983, ApJ, 266, 806 

Landolt, A.U., 1992, AJ, 104, 340 

Lenz, P & Breger, M., 2005, CoAst, 146, 53 

Levine, A. M., Brandt, N., Cui, W., Jernigan, J. G., Morgan, E. H., Remillard, R., Shirey, R. E. & Smith, D. A., 
1996, ApJ, 469, L33 

Lopes de Oliveira, R., Motch, C, Haberl, F, Negueruela, I. & Janot-Pacheco, E., 2006, A&A, 454, 265 
Lopes de Oliveira, R., Motch, C, Smith, M., Negueruela, I. & Ton-ejon, J.M., 2007, A&A, 474, 983 
Markert, T.H., Laird, FN., Clark, G.W., Heam, D.R., Sprott, G.F, Li, FK., Bradt, H.V., Lewin, W.H.G., Schnopper, 
H.W. & Winkler, RF, 1979, ApJS, 39, 573 



18 



G.E. Sarty et al. 



Mennickent, R.E., Pietrzynski, G., Gieren, W. & Szewczyk, O., A&A, 393, 887 

Menzies, J.W., Cousins, A.W.J., Banfield, R.M. & Laing, J.D., 1989, South African Astronomical Observatory 
Circular, 13, 1 

Motch, C, Belloni, T, Buckley, D., Gottwald, M., Hasinger, G., Pakull, M.W., Pietsch, W., Reinsch, K., Remillard, 

R.A., Schmitt, J.H.M.M., Trumpler, J. & Zimmermann, H.-U., 1991, A&A, 246, L24 
Motch, C, Haberl, R, Dennerl, K., Pakull, M. & Janot-Pacheco, E., 1997, A&A, 323, 853 
Okazaki, A. & Neguemela, I., 2001, A&A, 377, 161 

Press, W.H., Flannery, B.P., Teukolsky, S.A. & Vetterling, W.T., 1986, Numerical Recipies: The Art of Scientific 

Computing (Cambridge:Cambridge University Press) 
Pojmanski, G., 1997, Acta Astron, 47, 467 

Reichart, D., Nysewander, M., Moran, J., et al. 2005, II Nuovo Cimento C, 28, 767 

Reig, P, Fabregat, J. & Coe, M., 1997, A&A, 322, 193 

Reig, P & Roche, P, 1999, MNRAS, 306, 100 

Robinson, R.D., Smith, M.A. & Henry, G.W., 2002, ApJ, 575, 435 

Sarty, G.E., Kiss, L.L., Huziak, R., et al., 2009, MNRAS, 392, 1242 

Schmidtke, PC. & Cowley, A.P, 2006, AJ, 132, 919 

Smith, M.A. & Robinson, R.D., 1999, ApJ, 517, 866 

Smith, M.A. & Balona, L., 2006, ApJ, 640, 491 

Steele, I, Neguemela, I. & Clark, J., 1999, A&AS, 137, 147 

Stetson, PB., 1987, PASP, 99, 191 

Thackeray, A.D., Tritton, S.B. & Walker, E.N., 1973, Mem. R. Astron. Soc, 77, 199 
Torrejon, J.M. & Orr, A., 2001, A&A, 377, 148 

Tuohy, I.R., Buckley, D., Remillard, R.A., Bradt, H.V. & Schwartz, D.A., 1998, International Symposium on the 

Physics of Neutron Stars and Black Holes, 93 
Voges, W, 1992, 'The ROSAT all-sky X ray survey' in Environment Observation and Climate Modelling Through 

International Space Projects. Space Sciences with Particular Emphasis on High-Energy Astrophysics, 9 
Warwick, R.S., Marshall, N., Eraser, G.W., Watson, M.G., Lawrence, A., Page, C.G., Pounds, K.A., Ricketts, M.J., 

Sims, M.R. & Smith A., 1981, MNRAS, 197, 865 
Waters, L., Pols, O.R., Hogeveen, S.J., Cote, J. & van den Heuvel, E.RJ., 1989, 291, 775 
Wegner W., 1994, 270, 229 

Wen, L., Levine, A.M., Corbet, R.H.D. & Bradt, H.V., 2006, ApJSS, 163, 372 
Young, PJ., Roberts, W.H. & Sebo, K.M., 1999, ASPC, 172, 115 
Zorec, J., Fremat Y. & Cidale, L., 2005, A&A, 441, 235 



8.3 



> 




0.2 0.4 0.6 0.8 

HJD-2450000 



